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Abstract 

Cyclotron maser and plasma wave growth which results from electrons accelerated 
in magnetic loops are studied. The evolution of the accelerated electron distribution is 
determined by solving the kinetic equation including Coulomb collisions and magnetic 
convergence. We find that for modest values of the column depth of the loop the growth 
rates of instabilities are significantly reduced and that the reduction is much larger for 
the cyclotron modes than for the plasma wave modes. The large decrease in the growth 
rate with column depth suggests that solar coronal densities must be much lower than 
commonly accepted in order for the cyclotron maser to operate. The density depletion 
has to be similar to that which occurs during auroral kilometric radiation events in the 
magnet osphere . 

The resulting distributions are much more complicated than the idealized distribu- 
tions used in many theoretical studies, but the fastest growing mode can still simply be 
determined by the ratio of electron plasma to gyrofrequency, U = cj p /O e . However, the 
dominant modes are different than for the idealized situations with growth of the z-mode 
largest for U< 0.5, and second harmonic x-mode (s=2) or fundamental o-mode (s=l) the 
dominant modes for 0.5 < C/ <1. The electron distributions typically contain more than 
one inverted feature which could give rise to wave growth. We show that this can result in 
simultaneous amplification of more than one mode with each mode driven by a different 
feature and can be observed, for example, by differences in the rise times of the right and 
left circularly polarized components of the associated spike bursts. 

Subject earth: aurorae - hydromagnetics - particle acceleration - radiation mech- 

anisms - Sim: corona - Sun: radio radiation - wave motions 
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I. Introduction 

Rapidly fluctuating, short duration (1-100 ms) spike bursts with very high brightness 
temperatures (>10 12 K), high degrees of circular polarization (often up to 100%), and 
narrow bandwidth (Aw/u <0.01-0.1) have been observed at microwave frequencies from 
the solar corona (Droge 1977, Slottje 1978, Zhao and Jin 1982, Stahli and Magun 1986) and 
from the dwarf M star AD Leonis (Lang and Wilson 1986). Similar emissions, but of longer 
duration have been observed from the stellar system Algol (Lestrade et al. 1988) and other 
red dwarf (dMe) flare stars (Bastian and Bookbinder 1987, Kundu and Shevgaonkar 1988). 
The auroral kilometric radiation (AKR) in the Earth’s magnetosphere (e.g. Gurnett 1974, 
or more recently Calvert 1981), Jovian decameter radiation (DAM), and Saturnian (Carr, 
T.D., Desch, M.D., and Alexander, J.K. 1983) and Uranism (Gurnett et a/.1986) kilometer 
radiation also have similar characteristics. A common feature of all these objects is that 
they contain energetic electrons in magnetic loop structures and the frequencies of the 
bursts are on the order of the electron gyrofrequency for the magnetic field strengths and 
electron energies expected in their magnetospheres. 

The properties of the emission indicate that the radiation mechanism is not sponta- 
neous, but is coherent emission from electrons within the magnetic structure. Wu and Lee 
(1979) were first to recognize that the cyclotron resonance is important for nonrelativistic 
or midly relativistic electrons and that the dynamics of electrons of these energies in the 
Earth’s auroral zone can lead to an inverted momentum distribution which gives rise to 
amplified cyclotron emission (also known as cyclotron maser emission). This has become 
the accepted mechanism of AKR and has been suggested to explain solar and stellar spike 
bursts (Holman, Eichler, and Kundu 1980, Melrose and Dulk 1982a). Since then much 
theoretical work has been done to advance our basic understanding of cyclotron maser 
emission (Lee and Wu 1980, Hewitt, Melrose, and Ronnmark 1982, Sharma, Vlahos, and 
Papadopoulos 1982, Wu et al. 1982, Freund et al. 1983, Melrose, Hewitt, and Dulk 1984, 
Wu 1985). These works have established the relationship between the amplified cyclotron 
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emission and idealized inverted electron distributions (e.g. so called sin-N loss cone, hol- 
low beam, and DGH distributions) which are used to make the calculations tractable. It 
has been shown through these studies that, in all cases, the ratio of the electron plasma 
frequency to gyrofrequency, U = ui p /Q. e , determines the fastest growing mode. The basic 
results are that for small U < Uo < 1 the fundamental (s=l) x-mode is the fastest growing 
mode, for slightly larger U <1, the emission is fundamental z-mode or o-mode, and as U 
is increased further the dominant mode changes from z-mode to second harmonic (s=2) 
x-mode to second harmonic o-mode. For example, for the sin-N loss cone distribution 
Uo = 0.3 (see e.g. Melrose, Hewitt, and Dulk 1984, Fig. 1). We parenthetically add that 
the cyclotron maser mechanism has also been suggested to be the underlying physical 
process responsible for metric and decimetric continuum bursts from the sun (Winglee 
and Dulk 1986), localized heating of the corona (Melrose and Dulk 1982b), and particle 
acceleration (Sprangle and Vlahos 1983). 

In general, one requires a knowledge of the momentum distribution of electrons for 
determination of the characteristics of the instabilities and of the associated emission prop- 
erties. However, with the exception of the description of the electrons in the auroral zone by 
Wu et a/.(1982) or in coronal loops by White, Melrose, and Dulk (1986, hereafter WMD86), 
all of the theoretical discussions have assumed some idealized electron distribution with 
little attention to the details of how such inverted distributions are actually realized. The 
standard paradigm is that electrons are accelerated (or heated) in a magnetic loop with 
converging field structure. Then, after reflection their distribution is inverted due to an 
excess of perpendicular momentum ( df/dp± > 0) with the form of the inverted distribu- 
tion given by an idealized analytic form. The inversion gives rise to coherent emission, in 
either the x-mode, o-mode, z-mode, or some combination of these depending on the value 
of U in the source. These investigations have merely shown that this simple model can 
account for the general aspects of the observed spike bursts. 

There are, however, a few shortcomings of the theory in its simplest form. For example: 
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i) Typically, very high degrees (nearly 100%) of circularly polarized radiation are 
expected in the sense of the fastest growing mode with right (left) circularly polarized 
emission corresponding to the x-mode (o-mode), but for the majority of solar millisecond 
spike bursts the observed circular polarization is between 25 and 30% (Benz 1986). 

ii) The spike bursts generally occur during the rise and maximum phase of the impul- 
sive microwave emission (Slottje 1978, Zhao and Jin 1982, Zhao 1983) while formation of a 
loss cone distribution occurs after trapping and many bounce periods (see e.g. Aschwanden 
and Benz 1988). 

iii) The growth rates (amplification factors) computed using the distribution measured 
in aiiu by the S3-3 satellite in the Earth’s auroral zone (which shows both loss cone and 
hole structures) axe much lower than those found from the idealized distributions and are 
not large enough to account for the observed levels of AKR (Omidi and Gumett 1984). 
However, as pointed out by the authors in their article, the lower growth rates may be 
attributable to relaxation of the distribution before measurements were made or to the 
difficulty of measuring sharp discontinuities. 

All of these indicate that the actual distributions of electrons present in physical mag- 
netic loops are probably much more complicated than the idealized distributions assumed 
in the above studies. In order to evaluate the importance of departures from the idealized 
situations it is necessary to consider the dynamics of accelerated electrons in magnetic 
loops. WMD86 examined electron propagation effects in a model where electrons are 
impulsively heated at the apex of a converging magnetic loop, but did not include any 
stochastic scattering or energy loss processes e.g. Coulomb collisions. Furthermore, they 
compute the growth rates at times when the electrons are moving predominantly down- 
ward, while growth rates obtained after reflection are of primary interest because this is 
when the loss cone feature and maser action are expected to appear. 

The purpose of this paper is to remove these limitations of the model by considering 
the transport and evolution of the electron distribution in detail. In addition to mag- 
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netic mirroring, other physical interactions can be important. The relevant processes to 
consider are determined by the energy of the electrons and the properties of the ambient 
plasma. For example, energy loss and pitch angle scattering due to Coulomb interactions 
axe the dominant processes for electrons of energy 10 to 10 3 keV in flaring solar coronal 
loops (densities 10 9 to 10 n cm -3 and magnetic field strengths ~ 10 2 to 10 3 G). These and 
other interactions with the ambient plasma tend to smooth out the inverted features in the 
electron distribution and therefore reduce the growth rates and alter emission properties 
of the cyclotron maser. In this article, we determine the evolution of the energetic electron 
distribution including Coulomb interactions as well as magnetic mirroring and then exam- 
ine the subsequent effects on the growth rates of the cyclotron maser and plasma waves 
for different values of the field convergence and plasma density. We do not calculate the 
evolution of the wave spectrum and therefore do not include any nonlinear (or quasi-linear) 
effects on the electron distribution. 

In §11, we review the wave modes of interest and the relation between these modes 
and the momentum distribution of electrons and also discuss the method used to compute 
the growth rates. In §111, we give the kinetic equation which determines the evolution of 
electrons in magnetic loops and describe its solution. Then, in §IV, we use this solution to 
compute the growth rates for cyclotron maser and plasma wave emission for a wide range 
of physical conditions. We discuss our results in §V. 

II. Wave Properties and Growth Rate Formulae 
a) Electron Cyclotron Maser 

The electron cyclotron maser is associated with induced emission of the magnetoionic 
modes of plasma theory. The modes are the solutions to the dispersion equation for a cold 
plasma (collision frequency <C u; p , fl e ) in a uniform magnetic field at frequencies where the 
motion of the ions can be neglected (see e.g. Stix 1962, p.38). The wave properties are 
completely determined by the two ratios X = u> 2 /u> 2 and Y = Ct e /u>, where u> is the wave 
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frequency, and the angle 9 between the wave vector k and the direction of the magnetic 
field. There are two solutions to the dispersion equation at each frequency and we adopt 
the conventions of Melrose (1980b, p. 258) to describe these modes and their growth rates. 
The wave properties, namely the index of refraction n a = Kc/u, the longitudinal part of 
the polarization K „ , and the axial ratio of the polarization ellipse T<r, of mode a (-1 for 
x-mode, +1 for o-mode) are 


T„ — Y cos 6 
XY sin# TV 
Ka ~ \-X T„-Y cosO ’ 

T„ = -<r{x 2 + 1) 1/2 - x , 

Y sin 2 9 

X 2(1 — X) cos 9 

The resonant (n£ = oo) and cut-off ( n 2 a = 0) frequencies determine the frequency 
where these modes propagate. The waves propagate above their cut-off frequency 
is equal to the plasma frequency, u> p , for the o-mode, and a frequency 


( 1 ) 

( 2 ) 

(3) 

(4) 


range 

which 


^[l+(l+W 2 ) ,/2 ] . 


(5) 


for the x-mode. For frequencies below the cut-off, but larger than the resonant frequency 
the waves decay. The resonant frequencies are 

u>±(0) = ~ (l + £7 2 ± [(i + U 2 ) 2 — AU 2 cos 2 0] 1/2 ) , (6) 

where + is for the x-mode and — for the o-mode. Note that the usual sign convention, 
which we use here, is opposite to the usual sign convention used for a above. Below the 
resonant frequencies the modes change character. The x-mode becomes the z-mode which 
propagates in the frequency range u> x — Cl e < u < u>+. The z-mode can not escape a 
physical plasma because it will eventually be absorbed in a lower density region where its 
frequency equals the resonant frequency oj+. Below w_, the o-mode becomes the whistler 
which will not be of interest here. 
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The electron cyclotron maser action is possible when the electron momentum distribu- 
tion is such that amplification (negative damping) of these modes occurs. Near resonances 
one must include the effects of finite temperature of the plasma. This changes the disper- 
sion relation and wave properties given here for the cold plasma modes. Below, we discuss 
waves near the resonances further and their relevance to the maser problem. 


i) General Growth Rate Equations 

The growth rate r <T (k) of a particular wave (specified by mode, frequency, and angle) 
is given by (see e.g. Melrose 1980b, p. 275) 


oo - 

r *( k ) = ^2 / - s ^th - k \\P\\ c ) 


s£l e d d \ 


where the transition rate A* is 


^(k,p) = 


(2-Kecfi j.) 2 


K a sin 8 + (cos -n <r /?||)T (T 
sin 8 


Uz) + J',(z) 


+ T*)\ n.0 ± «u9 

J s is the bessel function of order s, its argument z = 'yn^^ sin0/F, and 


d{um a ) 


= 1 + 


XYT a cos 8 
2(7V — Y cos#) 2 


(1+JT)(1-T. 2 ) \ 

* (1 - X)(l + T2) ) ■ 


Including only the s = n term in equation (7) gives the growth rate of the n th harmonic of a 
particular mode a (e.g. a = —1, s = 1, yields the growth rate of the fundamental x-mode). 
The resonance. condition given by the 6 function in equation (7) reduces the integral to a 
one-dimensional integral along an ellipse in velocity space (hyperbola in momentum space) 
which is called the resonance ellipse. This condition severely limits the electron energies 
which contribute to the growth of a particular wave. Furthermore, the cyclotron maser is 
driven by nonrelativistic electrons therefore the argument of the bessel f un ction is small 
z 1 and J a (z) ~ z* so that the growth rates are largest for the first few harmonics. 
Therefore, the growth rate of a wave is dominated by a particular term (harmonic) in the 
sum so that it is useful to discuss the growth at different harmonics separately. 
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ii ) Equations for Nonthermal Electron Distributions 
Given the electron distribution /, calculation of the growth rates for the various modes 
is straightforward. As described in the next section, we will use the kinetic equation to 
obtain a solution for the electron distribution / as function of E and The transformation 
of variables necessary to apply the growth rate formula of equation (7) is obtained from 
f(E,n) = 7 2 /3m 3 c 3 / (p|| , p j_ ) , which gives 


rrc k) 


= S iEk n 


COS 9 \ 




( 10 ) 


M— MR 


f dE + r^ cose ~ fi) '^\ v - ^ 

where /j,r = (1 - sY/~f)l@n a cos 9. In general, / is the sum of the distribution of thermal 
background plasma and the nonthermal electrons whose distribution is obtained from 
equation (20) below. The background thermal plasma (with total density n 0 ) gives a 
negative contribution to the growth rate which we will consider separately, so that, unless 
otherwise specified, in what follows / stands for the normalized (f f dEd^.—\) nonthermal 
component of the distribution with total density n a . After further simplification, we find 
that the ratio of the growth rate to the gyrofrequency is 

Ei flEl f dE\[K 9 sm9 + {cos6 -n a fo R )T a )J a {z) + P ± n a sm0J' s {z) 

il e n 0 J I 




dE 0* 




ff(E 

'■0 J 


V t 2 /? 


( 11 ) 




( 12 ) 


B °( k) = "I™* *1 cos *l ((! + T D + * 

The calculation of the growth rate for any wave is straightforward using these equations 
if the distribution / is known. However, given / the determination the wave having 
the largest growth rate is not. In §IV below, we discuss an efficient method for this 
determination. 

We have tested this method of growth rate calculation and our expression for the 
growth rate given in equations (11) and (12) by computing the maximum growth rates for 
the sin-N loss cone distribution 


/<A«> = e*p(-0 2 /2#) { jta'W**.). « S « 


c 


(13) 
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where f(E, fi) = 7 3 /3/(/?, a), N = 6, a c = 30°, and fix = \/J?TJrnc* (with T = 10 8 K), and 
comparing our results to those of Melrose, Hewitt, and Dulk (1984, Fig. 1) and Aschwanden 
and Benz (1988, Fig. 5). We confirm these calculations. However, we note that our 
convention (see also Sharma, Vlahos, and Papadopoulos 1982) is to combine the factor of 
the ratio of the densities with the growth rate, while these previous works fix the density 
ratio (10 -2 ) for a particular value of U (0.1). There is then a difference of a density factor 
in the expressions for the growth rates, but the relative growth rates of the given modes 
are identical for any value of U . 


b) Plasma Waves 

In addition to maser emission, the distribution of accelerated electrons is expected 
to generate longitudinal waves (or plasma waves), especially when the electrons are pre- 
dominantly down going and therefore beam-like (see e.g. Hamilton and Petrosian 1987). 
In a magnetized plasma, these waves occur near a resonance of the magnetoionic modes 
with the value of U determining the resonance u>+(0) or u_(0). Near these frequencies, 
the wavelength becomes smaller than the Debye length and thermal corrections to the 
wave properties Eire then necessary (Melrose 1980, p. 265). For U < 1, the resonance 
at u = u>-(9) yields a longitudinal mode which is usually called the slow plasma wave. 
For U > 1, the longitudinal mode occurs at the resonance u) = u>+(0) and is called the 
generalized Langmuir wave. 

The growth rate of generalized Langmuir waves ( U > 1) was examined by McClements 
(1987) to determine the stability of the steady state solutions of Leach and Petrosian (1981) 
and in WMD86 the growth rate of the slow plasma wave ( U < 1) is evaluated for down 
going electrons. One finds in both of these papers that the maximum growth rate occurs 
at 9 = 0. Therefore, we evaluate the growth rate of plasma waves when the wave vector 
is along the magnetic field (9 = 0). This simplifies the analysis because, regardless of the 
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value of U, the frequency of these waves is given by 




= u)p + 3 k 


2 2 


( 14 ) 


where v\ = kT t jm t is the thermal velocity of electrons in the ambient plasma of temper- 
ature T e . For this angle of propagation only the s = 0 term in the sum of equation (7) is 
nonzero and the growth rate can be written 


Up n 0 n,ui J L 


+ 


7 (1 - n 2 ) d 


\\ 7 2 /? / 


( 15 ) 


M =i jmf3 


ft 2 ft dfi 

Thermal damping must also be considered for the plasma modes. This provides a negative 
contribution of 

■p nr, . i 

(16) 


- = -l/l-i «p(-l/2n?^) , 

ujp V 2 u >i n i 


so that the total growth rate is given by the sum of equations (15) and (16). 


III. Electron Dynamics and Distribution 

As discussed above, in most models of maser emission it is argued that the electron 
distribution is given by an idealized functional form such as the sin-N loss cone (see eq. 
[13]) and the emission properties are studied using growth rates computed using this ide- 
alized distribution. In this section, we present a more realistic model in which a specified 
electron distribution is injected in a magnetic loop. The electron distribution resulting 
from subsequent evolution and transport of the electrons in the magnetic loop which de- 
termines the growth rate of the modes is found by solving the electron kinetic equation 
using the Fokker-Planck formalism. This allows us to determine how the strength and 
degree of convergence of the magnetic field and the plasma density effect the maser action. 

a) Kinetic Equation 

The propagation and evolution of accelerated electrons in magnetic loops is deter- 
mined by the kinetic equation. For electrons of energy 10 keV to 1 MeV in an ambient 
plasma with density 10 9 to 10 14 cm -3 and magnetic field strength 10 2 to 10 3 G (conditions 
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believed to occur in the coronae of the sun and flare stars with ft e 0.3 to 3 GHz) the 
dominant processes are magnetic mirroring, Coulomb collisions, and wave-particle inter- 
actions. Observations allow us to estimate with some degree of confidence the expected 
effects of the first two processes, but we lack direct knowledge of the amount of plasma 
waves that are present in the flare plasma. The minimal level of plasma waves is that 
produced by the instabilities which we are about to evaluate. As mentioned above, in- 
clusion of the back reaction of the waves generated by electrons on their propagation is a 
nonlinear problem which increases the complexity of the situation beyond the scope of this 
paper. The steady state problem for U » 1 (only Langmuir waves) including the effects 
of collisions and the nonlinear effects of waves was treated in our earlier paper (Hamilton 
and Petrosian 1987). In this paper, we ignore the uncertain effects of the plasma waves 
so that the evolution of the electron distribution in a cold plasma (jfcT < typical electron 
energy) is determined by (see e.g. Hamilton, Lu, and Petrosian 1990) 


*L = -pcuV + 0 c d -hl A ( H zj£1 A + ±JL(t\ + _1_JL 

dt P P 0a P ds d/i V 2 V + A 0 dE\p) + A 0 /?V dfi 

(17) 


(1-V) 


d£ 

dfil 


+ S 


where A 0 = (10 24 cm)/no(s) In A, n<j(s) is the background plasma density, In A(« 20) is the 
Coulomb logarithm, and S(E , /*, s, t) is the source. Here we have assumed that the electron 
gyroradius is much smaller than other relevant length scales (i.e. A 0 or (d lnn 0 /<is) _1 ) so 
that the electron distribution function f(E,n,s,t ) can be described by four variables, the 
particle kinetic energy E (in units of me 2 ), pitch angle cosine /i, position along the magnetic 
field s, and the time t. 


b) The Source Term 

The existence of nonthermal electrons with energies of 10 to 100 keV in solar coro- 
nal loops is well established from studies of impulsive hard X-ray bursts (see e.g. Kane 
et a/. 1980). The electrons are characterized by a power-law energy spectrum f(E) oc E~ 6 
with 2.5 <6 i$7.5, but the angular distribution is not very well known (see e.g. McTiernan 
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and Petrosian 1990). Electrons of these energies (10 to 100 keV) accelerated on open field 
lines (or on closed field lines of very large radii) are also believed to be the source of type III 
radio bursts (Goldman and Smith 1986). It is therefore natural to assume that under cer- 
tain circumstances such accelerated electrons when reflected before reaching the transition 
region can be a source of free energy to drive the cyclotron maser. This assumption is 
supported by the association of the impulsive hard X-ray and spike emissions (see Benz 
and Kane 1986). Similar energetic particle spectra are expected to be present in other 
flare stars. 

Hence, we will assume a source function 

S(E, h ,s,t) = w ^j g -(Eo/E)‘eM-^/Ll)m . (18) 

This corresponds to the impulsive acceleration of electrons with an isotropic pitch angle 
distribution and power-law index <5 within an acceleration region of characteristic size L a . 
The distribution is normalized such that the density of nonthermal electrons of energy 
greater than Eq at s = 0 is n a . 

c ) The Ambient Plasma 

We assume, as is generally done, that the above distributions of electrons are injected 
into a symmetric coronal loop of half length L c with apex at 3 = 0. As is evident from 
equation (17), in addition to the source term 5, we also need to specify the plasma pa- 
rameters A 0 and d In £ /ds and their variation along the loop which are obtained from the 
variations of plasma density n 0 (s) and field strength £(s), respectively. For the purpose 
of solving equation (17) we only need the ratio £(s)/£<j where Bo is the field strength at 
s = 0 and not the absolute magnitude of the magnetic field. Following WMD86, we adopt 
the simple model 

B(s)/Bo = 1 + &* / Lrg (19) 

which corresponds to the variation of magnetic field strength in the far field of a dipole. 
We further ass um e that the density in the corona is constant and increases rapidly below 
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the transition region with a scale height much shorter than Lb or L c so that particles 
which reach there are quickly thermalized and do not contribute to the maser emission. 

The extent of the loss cone distribution is determined by the ratio Lb/L c or by what 
we call the mirror ratio r m = 1 + L\jL\ which is defined such that, in the absence of 
collisions, electrons injected (at s = 0) with pitch angles greater than a c = arcsinfl/^r^) 
are reflected above the transition region. On the other hand, the degree to which the 
collisions smooth out loss cones is determined by the column depth to the transition region 
Ntr ~ TIqL c . 

d) Solution of the Kinetic Equation 

We have solved the kinetic equation (17) using the numerical code described in 
Hamilton, Lu, and Petrosian (1990) for the impulsive source given by equation (18) with 
L a = 0.1£ c and <5 = 5. For the ambient plasma, we use r m = 1, 2, 5, and 10 and column 
depths to the transition region N tr = 0, 1.91 x 10 17 , 4.74 x 10 18 , 1.88 x 10 19 , 7.37 x 10 19 , 
and 4.36 x 10 2O cm -2 . Since the column depth an electron of energy E traverses before 
stopping is 

N(E)=n 0 X 0 E 2 /(E+l) . (20) 

the above column depths ( Nt r ) correspond to the stopping depths for electrons of initial 
energies E = 0,1, 5, 10, 20, and 50 keV, respectively. Strictly speaking, no plasma modes 
can be expected for Ntr = 0 or for zero density plasma. This case represents the idealized 
situation for which collisions can be neglected and provides a useful comparison for the 
finite density situations where collisions with the ambient plasma become increasingly 
important (cf. White, Melrose, and Dulk 1983). This case will also be useful for comparison 
with the results of WMD86 who ignored collisions. 

It should be noted that the results presented below are more general because of the 
simple scaling properties of the kinetic equation (17). The dynamics are completely deter- 
mined by the column depth N tr = n 0 L c and the mirror ratio r m . Reseeding the variables 
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in the kinetic equation according to the transformations V = rjL for all the lengths in 
the problem, n' = n/77 for the densities, and t' = r)t for the time, leaves the equation un- 
changed. Therefore, if the kinetic equation is solved for a given half loop length L c , density 
no, and mirror ratio r m , the solution for this mirror ratio and column depth N tr = n 0 L c 
is obtained for any half loop length L' c = t]L c and density n ' 0 = n 0 /rj by transforming the 
time t' = rjt, position s' = 77s, and L' a = tjL s . 

Some of the results from these calculations axe shown in Figures 1 and 2 where we 
present contour plots in velocity space (/?||, 0 ±) of the distribution function at four times 
(described below) and at two locations s = 0 and s = L c / 2 , respectively. The mirror ratio 
r m = 5 corresponds to L c /Lb = 2 or critical pitch angle a c = 26° (/x^ = 0.8). On each 
figure, the left, middle, and right panels axe for column depths corresponding to stopping 
depths for electron energies E = 0, 10, and 20 keV, respectively. In the absence of collisions, 
electrons of initial energy E, (or velocity /?,) and pitch angle cosine /x, will be reflected 
after a time Tb = 7tTb/2^,c>/ 1 — /x?. Since we consider electrons of energy E > E 0 (see 
eq. [18]) and because only electrons with m < \i c (or a< > a c ) are reflected, the maximum 
value of interest for Tb is Tg az = 7tTb v /r^/2/3 0 c. We use a minimum electron energy 
Eq = 10 keV, and display from top to bottom on each figure the distribution at times 
*1 = ylfl ' BX , <2 = \T$ ax , t 3 = T5* ai , and U = f Tg ax , respectively. 

Vertical and horizontal cuts on this figure show where df/d |/?||| > 0 and df/d 0 ± > 0 
which axe the features of the electron distribution that give rise to amplified emission. For 
example, in Figure 1 for the model E = 0 along the line /3|| = 0.3 at time <2 there is one 
region and at <3 there are three regions where df /d() j_ > 0. Regions of df/d( 3 ±_ > 0 are 
also seen in Figure 2 at these times. Positive gradients in parallel velocity (d//d|/?||| > 0) 
axe most clearly seen in Figure 2 at , for all values of the column depth. 

We have shown contour plots in velocity space even though we calculate the distribu- 
tion f(E, /j.) in energy and pitch single cosine space and the growth rate formulae are given 
in terms of momentum space variables (see eq. [7]). We do this because the phase space 
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defined by (/?||,/?_l), unlike the other variable choices, is finite and because it is common 
practice to display these contours in velocity space. Notice that for the energies involved 
there will be very little difference between the shapes of these contours and those made in 
momentum space so that features found in the velocity distribution here axe also present 
in the momentum distribution. No contours are shown for f3 < 0.2 because Eq = 10 keV 
is the lowest energy considered in our calculations. 

Comparison of the two figures shows that the distributions at the top of the loop are 
markedly different than distributions midway down the loop. This would indicate that the 
emission properties such as polarization, frequency, angle, and brightness temperature will 
vary along the loop. Comparison of the contours for different column depths shows that 
the effect of the collisions is to smooth out the sharp features of the distribution which are 
the source of instabilities so that even for modest column depths (lV tr >7.37 x 10 19 cm -2 ), 
the unstable features of the distribution are only seen on the lowest contours and at early 
times. 

Note that these more realistic distributions are far from the idealized assumed dis- 
tributions mentioned earlier. However, parts of these distributions resemble parts of the 
idealized distributions. For example, the distributions in Figure 1 at have loss cone 
features which are qualitatively similar to the sin-N distribution, and the distribution at 
t 3 for N tr = 0 in Figure 2 shows isolated contours centered at (/?|j rsj 0. 1 , Al 0.3) which 
resemble the so called hollow beam distribution. 

IV. Growth Rates and Analysis 

As mentioned above, when determining the dynamical evolution of the distribu- 
tion we need to specify the rate of convergence (or divergence) of the magnetic field, 
dlnB/ds , but not its magnitude. However, as stated in §11, the most important pa- 
rameter in determining the growth rates is the ratio of the plasma to gyrofrequency 
U = (n/10 1 °cm _3 ) 1 / 2 (320G/B). Therefore, to completely specify a model, in addition 
to the mirror ratio and the column depth, we must specify the value of U at the top of 
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the loop. We should also consider the effects of finite temperature even though these are 
unimportant for determining the particle distribution as long as kT <C •E'o • The primary 
effect here is thermal damping of the waves by the ambient plasma. However, for magne- 
toionic waves, damping is unimportant provided that the above inequality is valid so that 
particle velocities fax exceed the thermal velocity v e . 

a) Cyclotron Maser 

We need to compute the growth rates of x-mode, o-mode, and z-mode waves for a 
particular distribution /. In order to sample the frequency space (u,0) efficiently, we 
use the method outlined in Aschwanden (1989). For a given mode cr and ( u>,6 ) pair, the 
resonance condition represented by the 6-function in equation (7) 

u> — — k\\0\\c = 0 (21) 


defines an ellipse in velocity space (0\\,0±), 

0*11 - 0 o ) 2 


0 2 . 

+ — = 1 


F 2 (l — e 2 ) V 2 


( 22 ) 


with the center (J3 0 , 0), semi-major axis V, and eccentricity e given by 

e 2 

*° “ n a cos 6 ’ 

V = [l-Ce/n^costf) 2 ] 1 / 2 , ( 23 ) 

e = [1 + {sY(n a cos#) 2 ] -1 / 2 . 

Conversely, each point in velocity space ( 0 ||,/?_l) can be identified with the top of a reso- 
nance ellipse with 0q — V = f3±, and frequency u = sn e 7 (l — 0\). The last relation 
can be obtained from equation (21) using the expression fc|| = k cos 9 = cos d/c and 
the fact that at the top of the resonance ellipse e 2 = 0\\n„ cos 6 and 1 — e 2 = 7 _2 /(l — /3 2 ). 
The frequency space grid is set-up by first finding the frequency range defined by the ve- 
locity space grid and then subdividing this range. The range of angles for each frequency 
is determined by the requirement that the resonant condition be satisfied (eq. [21]). 
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Using the formulae of §11, results of §111, and this method of calculation, we evaluate 
the growth rates of the x-mode (s = 1,2), o-mode (s = 1,2), and z-mode (s=l) for three 
values U t0 p = 0.1, 0.5, and 1.0. The results of our calculations are shown in Tables 1 (top 
of the loop) and 2 (at 3 = L c /2) at times t 2 = Tg ax and t 4 = § Tg ax for three values 
of the mirror ratio r m . Note that the growth rates are functions of frequency and angle 
and that we list the maximum values of the growth rates of each mode here. The mode 
having the largest growth rate for a particular time and model is bold faced. We draw the 
following conclusions from the results listed in these tables. 

1. As the column depth of the plasma in the loop increases the growth rates decrease. 
This is a general trend which is independent of the values of r m or U. Even for a modest 
column depth of 1.88x 10 19 cm -2 the growth rates are reduced by 2 to 4 orders of magnitude. 

2. The mode with the largest growth rate does not remain the same throughout the 
entire loop because for converging magnetic field geometries the ratio U and the electron 
distribution vary along the loop. 

3. Although the distributions are complicated and vary along the loop, there are some 
general trends in the results. These trends are in general different than those obtained from 
assumed distributions or idealized situations. The fastest growing modes are z-mode for 
U <0.5, and second harmonic x-mode or fundamental o-mode for 0.5 < U < 1. Fundamental 
x-mode growth is not important, even for U < 0.3. The fastest growing mode also depends 
on the value of the mirror ratio r m . For example, for U top = 0.5, the mode having the 
largest growth rate at the top of the loop is the second harmonic x-mode when r m = 2 but 
changes to the fundamental o-mode for r m > 5. 

4. An important consequence of the complexity of the distribution is the presence of 
more than one inverted feature which occurs even at larger column depths but with reduced 
gradients. This leads to the possibility that each of these inverted features might drive 
a different instability. Since relaxation of one part of the distribution does not directly 
affect the other inverted features, the distribution can support growth in more than one 
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mode simultaneously. To illustrate this, we show in Figure 3 a contour plot of the electron 
distribution for the model with N tr = 188 x 10 19 cm -2 and r m = 5 at s = L c j 2 for time 
t 4 = ^Tg ax . Overlying the contours are the resonance ellipses for the o-mode (s=l), 
x-mode (s=2), and z-mode (s=l) corresponding to the frequencies and angles for these 
modes that produced the largest growth rates given in Table 2. From the table we see that 
T" > Tj > TJ, but as evident from Figure 3, these modes are driven by different features 
of the inverted electron distribution. Therefore, even though T® ~ 5rf, saturation of the 
o-mode (if it occurs) will not inhibit amplification of the x-mode. This is an important 
result since this provides a mechanism to explain saturated maser emission which is not 
very highly (~100%) circularly polarized. If the burst is resolved, a direct observational 
consequence is that the right and left circular polarizations will show different rise times 
due to the differences in the growth rates. It is clear that a more detailed analysis including 
the dynamics of the wave growth is required in order to quantitatively determine the time 
profile of the resulting spike bursts, but the qualitative behavior should agree with this 
description. 

b) Plasma Waves 

Because damping by thermal electrons is important, the growth rate does not simply 
scale with the ratio of the density of nonthermal to thermal electrons. The dispersion 
relation as well as the thermal damping rate are sensitive to the temperature of the ambient 
plasma. We compute the growth rate of plasma waves for n a /no = 10 -2 and 10 -4 for 
T = 10 6 and 10 7 K. In Table 3, we give the maximum growth rates of waves with wave 
vectors along the loop for times = $Tg ax and <2 = \Tg ax . As expected, increasing the 
column depth, increasing the temperature of the thermal plasma, or decreasing the ratio 
n a /n 0 reduces the growth rate. Notice that increasing the column depth does not reduce 
the production of plasma waves by as large a factor as was seen for the magnetoionic 
modes. Also note that the growth rates are largest at early times when the electrons are 
down going and the distribution has beam-like nonthermal features. 
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Relaxation of the distribution due to the generation of plasma waves results in the 
reduction of the column depth traversed by the forward moving electrons (see e.g. Hamilton 
and Petrosian 1987), but is not expected to influence the maser action since these processes 
involve different electrons (cf. WMD86). However, we compute these growth rates to 
emphasize the fact that whenever maser action occurs in this type of a physical system it 
is accompanied by the production of plasma waves. Furthermore, the plasma waves are 
generated at earlier times because they are driven by beam features which appear in the 
distribution before reflection of electrons occurs. However, plasma waves could produce 
density fluctuations which would modulate the maser emission. 

V. Discussion 

We have examined the cyclotron and plasma wave growth induced by electrons accel- 
erated in magnetic loops. We considered an initial accelerated electron distribution which 
is consistent with nonthermal electron distributions that produce the impulsive hard X-ray 
spectrum found in solar flares and is expected to be present in stellar flares as well. The 
evolution of the electron distribution was determined by solving the kinetic equation in- 
cluding magnetic convergence and Coulomb collisions with the ambient plasma. In general, 
the fastest growing mode is determined by the value of U, as for the idealized distributions, 
with z-mode dominant for U <0.5, and second harmonic x-mode or fundamental o-mode 
for 0.5 < U < 1. Growth of the fundamental x-mode is not significant which is in contrast to 
the idealized situation. The resulting distributions have a complex structure and typically 
have more than one inverted nonthermal feature. We showed that these features can drive 
different modes simultaneously and that this provides a mechanism to explain partially 
polarized emission. We note that it might also be possible for multiple features to give rise 
to simultaneous amplification of different harmonics of the same polarization (e.g. x-mode 
at s=l and s=2). Plasma wave growth is also shown to occur so that if cyclotron maser 
emission produces the spike bursts the density fluctuations which result from the plasma 
waves causes variations in U which could then modulate the emission. 
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When collisions are included maser action still occurs, but the growth rates are reduced 
significantly. The growth of plasma waves is also reduced, but by a much smaller amount. 
Under modest conditions of column depth N tr = 1.88 x 10 19 cm -2 the growth rates are 
typically reduced by 2 to 4 orders of magnitude for the magnetoionic modes and less than a 
factor of ten for the plasma waves. For nominal physical conditions ( L c ~ 10 9 -10 lo cm and 
no ~ 10 9 -10 u cm -3 ) in flaring solar coronal loops the column depth is this large or larger. 
Therefore, maser emission in such loops is expected to be severly limited. This is also 
true in larger ( L c > 10 lo cm) low density (n 0 < 10 9 cm -3 ) loops which extend into the upper 
corona because our results depend primarily on the column depth and to a lesser extent 
on the the ratio of magnetic field scale height to loop length. Under these conditions, 
this indicates that the spikey, rapidly varying bursts might be due to plasma emission as 
suggested by Kuijpers, van der Post, and Slottje (1981) and not cyclotron maser emission. 
However, the spike bursts should then be similar to type III bursts, but the characteristic 
signatures of the spike and type III bursts are very different (see Benz 1986 Fig. 2). In 
addition, density depletions are observed to occur simultaneously with the AKR (Benson, 
Calvert, and Klumpar 1980). This may indicate that similar processes are occurring in 
planetary, solar, and stellar magnetospheres. Our results suggest that in a solar flare 
region, where there exists a large range of magnetic field strengths, temperatures, and 
densities, maser action is possible only in the regions of lower density (and presumably 
higher magnetic field or temperature for pressure balance) and should occur primarily 
at the flare onset before chromospheric evaporation increases the density suppressing the 
maser action. This behavior explains why the spike bursts are usually seen only during 
the rise and maximum phases of the microwave emission during solar flares. 
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Table 1 


l0g (%^?) at t0p *°°P 
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Table 2 


lo '(% f ^) a,W2 


t: 

N c : 0 

2 

10 20 

h 

5 

0 10 20 

10 

0 10 20 

2 

0 10 20 

<4 

5 

0 10 20 0 

10 

10 20 

U t op — 

0.1 






U : 

0.08 

0.05 

0.03 

0.08 

0.05 

0.03 

XI: -6.35 

-6.46 -7.61 

-5.44 -5.74 -6.74 

-5.52 -5.46 -6.80 

-3.89 -7.24 -7.94 

-3.29 -6.14 -7.71 -3.52 

-6.44 -8.01 

X2: -3.62 

-5.17 -7.08 

-3.18 -5.10 -7.16 

-3.72 -5.54 -7.61 

-2.89 -6.31 -7.99 

-3.21 -6.84 -8.63 -3.86 

-7.42 -9.17 

Ol: -4.82 

-5.23 -6.83 

-4.03 -5.05 -6.69 

-4.53 -5.50 -7.02 

-2.81 -6.30 -7.77 

-3.32 -6.25 -7.85 -4.05 

-6.67 -8.15 

02: -6.22 

-6.36 -7.98 

-5.41 -6.53 -8.10 

-6.25 -6.98 -8.70 

-4.61 -7.50 -8.89 

-5.23 -7.93 -9.39 -5.88 

-8.59 -9.92 

Z 1 : - 2.53 

- 4.10 - 6.12 

- 1.88 - 3.71 - 5.75 

- 2.25 - 4.14 - 5.98 

- 1.53 - 5.25 - 6.93 

- 1.91 - 5.17 - 6.92 - 2.33 

-5.55 -7.21 

Utop — 

0.5 






U: 

0.40 

0.25 

0.15 

0.40 

0.25 

0.15 

XI: * 

* -7.54 

-4.67 -4.70 -5.54 

-4.42 -4.44 -5.57 

-9.49 -7.31 -8.54 

-2.96 -5.06 -6.75 -3.00 

-5.26 -6.86 

X2: -2.46 

-4.34 -5.54 

-1.94 -4.36 -5.53 

-2.03 -4.60 -6.12 

-1.37 -4.99 -6.68 

-1.91 -5.47 -7.11 -2.25 

-6.08 -7.89 

01: -3.51 

-4.10 -5.30 

-2.58 -3.69 -5.32 

-3.14 -4.02 -5.62 

-1.50 -4.97 -6.35 

-2.00 -4.85 -6.76 -2.66 

-5.27 -6.84 

02: -4.85 

-5.01 -6.42 

-4.00 -5.15 -6.72 

-4.86 -5.56 -7.31 

-3.25 -6.02 -7.51 

-3.85 -6.56 -7.93 -4.49 

-7.20 -8.52 

Z 1 : - 2.36 

- 4.07 - 5.56 

- 1.63 - 3.62 - 5.08 

- 1.73 - 3.79 - 5.18 

-1.48 - 4.84 - 6.35 

- 1.67 - 4.61 - 6.23 - 1.82 

-4.84 -6.39 

Utop = — 

1.0 

- 





U: 

0.8 

0.5 

0.3 

0.8 

0.5 

0.3 

XI: * 

* * 

-9.16 -9.09 -8.63 

-4.97 -4.86 -5.21 

* * * 

-8.92 -8.87 -8.70 -4.80 

-5.19 -6.53 

X2: - 1.33 

-3.55 - 5.31 

- 1.62 - 3.37 -5.17 

-1.87 -3.85 -5.59 

- 0.56 - 4.40 - 5.81 

- 1.24 -5.06 -6.45 -2.00 

-5.41 -7.14 

Ol: -2.12 

-4.21 -5.45 

-2.44 -3.53 - 4.53 

-2.58 -3.30 -5.11 

-1.10 -4.43 -6.27 

-1.29 - 4.32 - 5.61 -2.09 

-4.79 -6.53 

02: -4.37 

-4.76 -6.00 

-3.39 -4.61 -6.21 

-4.27 -4.91 -6.73 

-2.78 -5.82 -6.88 

-3.21 -5.96 -7.57 -3.92 

-6.63 -7.93 

Zl: -1.99 

-3.39 -6.13 

-1.86 -3.76 -5.00 

- 1.78 -3.82 - 4.89 

-1.85 -3.60 -6.34 

-1.77 -4.84 -6.16 -1.80 

-5.24 -5.89 


Values of the growth rate smaller than 10 -10 are represented by an asterick. 
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Table 3 





I'm- 


1 



2 



5 



10 



N c : 

0 

10 

20 

0 

10 

20 

0 

10 

20 

0 

10 

20 

t = 

ti 













n a /no 

T 













10~ 2 

10 6 

-3.22 

-3.61 

-4.41 

-2.03 

-2.76 

-3.78 

-2.81 

-3.32 

-3.88 

-3.57 

-3.94 

-4.39 

10- 2 

10 7 

-5.26 

-5.26 

-5.29 

-2.22 

i 

CO 

b 

00 

-4.01 

-3.54 

-3.70 

-4.06 

-5.24 

-5.28 

-5.34 

O 

1 

10 6 

-5.35 

-5.92 

-6.41 

-4.03 

-4.76 

-5.78 

-4.81 

-5.32 

i 

Oi 

bo 

00 

-5.69 

-5.94 

-6.40 

10" 4 

10 7 

-7.26 

-7.26 

-7.29 

-5.07 

-5.55 

-6.22 

-6.15 

-6.27 

-6.46 

-7.27 

-7.32 

-7.38 

t = 

*2 













n a /n o 

T 













10"2 

10 6 

-2.32 

-2.88 

-3.42 

-2.17 

-4.48 

* 

-2.54 

-3.82 

-5.75 

-2.85 

-3.92 

-5.66 

10-2 

10 7 

-2.84 

-3.22 

-3.66 

* 

* 

* 

* 

* 

* 

* 

* 

* 

10“ 4 

10 6 

-4.32 

-4.88 

-5.42 

-4.17 

-7.57 

* 

-4.54 

-5.82 

-7.96 

-4.85 

-5.92 

-7.66 

10" 4 

10 7 

-5.18 

-5.45 

i 

OI 

CO 

00 

* 

* 

* 

* 

* 

* 

* 

* 

* 


Values of the growth rate smaller than 10 -10 are represented by an asterick. 
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Figure Captions 

Fig. 1. Contour plots of the electron distribution function in the /3|| plane at the top of 
the loop for a mirror ratio of r m = 5. The column depths from left to right are N tr = 
0, 1.88 x 10 19 , and 7.37 x 10 19 cm -2 corresponding to stopping depths for electrons 
of energies E = 0, 10, and 20 keV, respectively. The times from top to bottom are 
lrpmax^ j^max , anc i | Tg ax . The contour levels are spaced by a factor of 10 

and axe the same for each plot. 

Fig. 2. Same as Figure 1, but at s = LJ2. 

Fig. 3. Contour plot of the electron distribution shown in the middle of Figure 2 with 
N tr = 1-88 x 10 19 cm -2 at time | Tg ax . Contour levels are now spaced by a factor 
of 5. The darker lines represent the resonant ellipses of the o-mode (s=l), x-mode 
(s=2), and z-mode waves at frequencies and angles which maximize their growth 


rates. 
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